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We study the spectral properties of a family of quantum antiferromagnets on two-dimensional 
(2D) lattices. This family of models is obtained by a deformation of the well-studied 2D quantum 
antiferromagnetic model of Affleck, Kennedy, Lieb and Tasaki (AKLT); they are described by two- 
body, frustration-free Hamiltonians on a three-colourable lattice of spins. Although the existence of 
a spectral gap in the 2D AKLT model remains an open question, we rigorously prove the existence 
of a gap for a subset of this family of quantum antiferromagnets. Along with providing new progress 
for the gap problem in AKLT-type antiferromagnets in 2D, this result has implications for the theory 
of quantum computation as it provides a family of two-body Hamiltonians for which the ground 
state is a resource for universal quantum computation and for which a spectral gap is proven to 
exist. 


I. INTRODUCTION 

Quantum antiferromagnetic spin lattice models pos¬ 
sess a rich structure, and exhibit many properties of 
quantum phases that lie beyond the standard Landau 
classification by symmetry breaking. In one-dimensional 
chains, quantum antiferromagnetic models were the sub¬ 
ject of Haldane’s conjecture: that there is a gap in 
the energy spectrum of integer-spin antiferromagnetic 
chains whereas half-odd-integer-spin chains are gapless. 
The analytically-solvable model of Affleck, Kennedy, Lieb 
and Tasaki (AKLT) for a one-dimensional spin-1 chain 
possesses a spectral gap, a ground state degeneracy 
determined by the boundary conditions, and a nonlo¬ 
cal ‘string-like’ order parameter CHS]; this model now 
serves as the canonical example of a one-dimensional 
spin chain possessing both symmetry-protected topolog¬ 
ical order [1] and a valence bond solid description. For 
higher-dimensional lattices, quantum antiferromagnetic 
models are expected to possess a similarly rich structure, 
but classifying these structures is a challenge, including 
even the rigorous identification of a spectral gap in nat¬ 
ural generalisations of the AKLT state to 2D lattices. 

Methods from the theory of quantum computation pro¬ 
vide a new approach to identifying and classifying the 
quantum order present in antiferromagnetic spin lattice 
models [5]. In particular, it has been shown for ID 
antiferromagnetic spin chains that their unique type of 
quantum order can be identified with an infinite local- 
izable entanglement length EH]. For the 2D AKLT 
model and variants on a number of spin-3/2 lattices 
and the spin-2 square lattice, an even stronger result 
has been proven: that the ground state of such mod¬ 
els possesses an entanglement structure that allows for 
universal quantum computation through local measure¬ 
ments alone [iDHnj. Quantum computation that pro¬ 
ceeds by single-particle measurements of a many-body 
entangled quantum state is called measurement-based 
quantum computation (MBQC). The universality of the 
AKLT model for measurement-based quantum compu¬ 


tation was proven by demonstrating that its ground 
states, on certain lattices, can be converted through lo¬ 
cal operations into a class of entangled states known 
as graph states —states that form the central resource 
for the theory of measurement-based quantum computa¬ 
tion [HIITS] . From the perspective of quantum comput¬ 
ing, there is keen interest in identifying and classifying 
the types of quantum order in such spin lattices that en¬ 
able quantum computation in this sense [niiniiiiiiiniiia- 
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In this paper, we show that methods from quantum in¬ 
formation theory, in particular the study of graph states 
and AKLT states for quantum computation, can be used 
in proving spectral properties of 2D antiferromagnetic 
spin models. In particular, we consider a family of 2D an¬ 
tiferromagnetic spin models obtained by a one-parameter 
deformation of the 2D AKLT model. This model was in¬ 
troduced in a previous work |24j . wherein it was shown 
using numerical methods that there is a region in pa¬ 
rameter space where the ground state of such models is 
universal for quantum computation. The deformation to 
the AKLT model we consider is a modification of that of 
Ref. In contrast to our case, however, the ground 

state of the model in Ref. [5S] tends to a Neel state (rather 
than a graph state) as the deformation becomes large. 

Here, we prove several rigorous results concerning the 
spectral properties of our family of antiferromagnets, 
with a focus on the key question of whether there ex¬ 
ists a spectral gap separating the ground state from the 
first excited state in the thermodynamic limit. While 
the AKLT model on which the model is based is not cur¬ 
rently known to be gapped, we prove that a subset of 
this family of models (not including the AKLT point) is 
gapped provided the ground state is sufficiently close to 
a graph state. Our work complements other studies that 
have numerically (rather than analytically) investigated 
the gap in other families of antiferromagnets based on 
the AKLT model [MIET]. 

One implication of our results is that this family of 
models has many of the desired properties for quantum 
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computation, in particular that the Hamiltonian is two- 
body and that there exists a finite region in the parameter 
space that is both computationally universal and gapped. 
Previous works have presented 2D models, based on the 
AKLT model, that are both gapped and universal for 
MBQC [HlEO]. However, these models do not possess 
the full 2D structure of the family of models presented 
here. The model in m is unitarily equivalent to a prod¬ 
uct of ID chains, and the model in [20] can be trans¬ 
formed via a unitary into decoupled regions of finite size. 

The paper is structured as follows. In Sec.|Hj we define 
the quantum spin lattice models that we investigate: the 
AKLT model, the graph state model, and our family of 
2D quantum antiferromagnets. We then investigate the 
spectral properties of these models in Sec. |HI[ including 
a rigorous proof of the existence of a spectral gap for a 
subset of these 2D quantum antiferromagnets. We sum¬ 
marise our results in Sec. US Several technical results are 
relegated to appendixes. 


II. AKLT MODELS, GRAPH STATE MODELS, 
AND A FAMILY OF 2D ANTIFERROMAGNETS 

In this section we define AKLT states and graph states 
in terms of their parent Hamiltonians. We note that 
AKLT states and graph states are special in that both 
are ground states of frustration free Hamiltonians (i.e., 
both AKLT states and graph states are the simultaneous 
ground states of all individual terms in their respective 
Hamiltonians). In addition, they are describable as pro¬ 
jected entangled pair states (PEPS) [28l and can be cre¬ 
ated from product states using constant-depth unitary 
circuits. 

We then define a class of 2D antiferromagnet mod¬ 
els obtained as deformed versions of the AKLT model, 
generalising the spin-3/2 models originally presented in 
Ref. [Ml- These deformed AKLT models have ground 
states that interpolate between an AKLT state and a 
graph state. 


A. Parent Hamiltonians for AKLT and graph 
states 

Here, we review parent Hamiltonians for graph states 
and AKLT states. Alternative equivalent definitions of 
these states in terms of tensor network states are detailed 
in Appendix [A] 

A parent Hamiltonian for a state \ip) of TV particles is 
a local Hamiltonian H = hi that has |^/>) as its unique 
lowest energy eigenstate. By local, we mean that each in¬ 
teraction term hi acts non-trivially on at most k particles, 
where fc is a constant independent of N. The Hamilto¬ 
nians that we consider in this paper are also frustration- 
free, meaning that the ground state of H is the lowest 
eigenstate of each interaction term hi individually. The 
minimum energy of each hi term may be set to zero by 


adding a multiple of the identity, and by replacing each hi 
with the projector onto its image we will not change the 
ground space of the Hamiltonian nor any crucial spectral 
properties such as the existence of a spectral gap. There¬ 
fore, for simplicity, in this paper we need only consider 
frustration-free Hamiltonians that are sums of projectors. 

1. AKLT states 

The standard AKLT parent Hamiltonian may be de¬ 
fined on an arbitrary graph G = (V, E) as follows. At 
each vertex v of the graph place a spin-S't, particle, where 
Sv = {dv — l)/2 and dy is the degree of vertex v £ V. 
The AKLT Hamiltonian is an antiferromagnetic interac¬ 
tion Hamiltonian given by 

jjA _ p,Stot=S.+S, ^ 

(*d> 

where Si is the spin of particle i, and is the pro¬ 

jection onto the spin-S' irreducible representation of par¬ 
ticles i and j. Defined in this way, has a unique zero- 
energy eigenstate, which we call the AKLT state |AKLT). 
Furthermore, the AKLT state is a zero eigenstate of each 
term in Q , implying that the Hamiltonian is frustration- 
free. We note that the terms in this Hamiltonian do not 
commute with each other, and so establishing the exis¬ 
tence of a spectral gap for such models is highly non¬ 
trivial. For two-dimensional lattices, it remains an open 
question in general whether this AKLT Hamiltonian pos¬ 
sesses a spectral gap in the thermodynamic limit. 

2. Graph states 

A graph state, and its corresponding parent Hamilto¬ 
nian, may also be defined on an arbitrary graph G = 
{V,E). In contrast to AKLT states, however, every par¬ 
ticle in a graph state is taken to have spin-1/2. (While 
we do not consider it here, it is possible to generalise 
to higher spin [IH]-) Following the original definitions 
[HI US], the graph state |G) corresponding to a graph G 
is defined as the output of a simple quantum circuit 

IG) = ((g)i+)J , (2) 

\eeE ) \vev / 

i.e., a qubit in the state |-|-) := ^(|0) -I- |1)) is 
placed at each vertex, and a controlled-Z gate CZ := 
exp(*7r|Il)(ll|) is applied between every pair of vertices 
connected by an edge. 

In contrast to AKLT states, graph states possess a par¬ 
ticularly simple parent Hamiltonian for which the spec¬ 
tral properties are much simpler. The graph state is the 
unique ground state of the Hamiltonian 

- K^), 


(3) 
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where Ki := stabilizer generators, Xi 

and Zi are Pauli X and Z matrices respectively acting 
on particle i, and where the product is taken over all 
neighbours rii of particle i. Note that the interaction 
terms {|(/ — Ki) : i = 1,..., N} pairwise commute. So, 
unlike the AKLT Hamiltonian, this cluster Hamiltonian 
is trivially diagonalisable and possesses a unit spectral 
gap. 

3. Relationship between AKLT and graph states 

Motivated by physical considerations, one generally 
views the Hamiltonian as being more realistic than 
PfC, as it involves only two-body interaction terms. The 
cluster Hamiltonian involves many-body interac¬ 
tions. Even for a one-dimensional chain (degree-two), 
the interaction terms Ki of each act non-trivially 
on three particles, and the bodiness of the interaction 
terms increases with the degree of the graph. Despite 
these somewhat artificial features, the cluster Hamilto¬ 
nian has commuting terms and a unit gap, and its 
simplicity makes it a useful tool for investigating spectral 
properties of more general models. 

Let G be three-colourable graph and let |AKLT) be the 
AKLT state defined on this graph. AKLT states on three- 
colourable graphs have the desirable property that they 
can be transformed into graph states on the same graph 
using a set of single-particle projectors m- We use this 
property to construct a continuous family of Hamiltoni¬ 
ans and their associated ground states that interpolate 
between these AKLT states and the corresponding graph 
states. 

Define three operators 

p-:=\S),{S\ + \-S)A-S\, (4) 

for c G {x,y,z} where |±S')c are spin-S” states satisfying 
5'^|±S')c = ±S'|±S')c and jS'^ is the single-spin operator 
for a spin-S' particle along the c axis. In other words, the 
operators are projections onto the two-dimensional 
subspace spanned by |±5')c. 

Now, let {ci G {x,y,z} : i = l...A^} be a three- 
colouring of G, i.e., Cj ^ Ci if * and j are neighbours. 
For each particle i in our system, we define with re¬ 
spect to Eq. 0. with S taken to be the spin of particle 
i. We will call the image of P^' the logical space of 
particle i and we will associate the physical spin states 
\S)c and I — S)c with the logical qubit states |0) and |1) 
respectively. Then the projectors P^' map the AKLT 
state to a graph state in the sense that 

N 

(g)i^“-|AKLT)cx|G), (5) 

i=l 

where |G) is a graph state on the graph G, where qubit 
i is encoded into the logical subspace as above. This 
result was originally derived in cni, and we provide an 
alternative proof using tensor networks in Appendix [B} 


B. A family of two-body antiferromagnets 

In this section we will construct a family of deformed 
AKLT-like antiferromagnet Hamiltonians, originally de¬ 
scribed in Ref. [ 21 ] for the special case of spin-3/2 models, 
that interpolate between the physically-motivated but 
more complicated AKLT Hamiltonian and a sim¬ 
pler model with a structure based on graph states. This 
Hamiltonian has the appealing property that it has only 
two-body interactions. 

For each site i, we define a one-parameter family of 
single-particle operators Di{S) acting on site i as 

D,{6) := {1 - 6)P,^'+ 61, (6) 

where 0 < 5 < 1. This operator satisfies Di{l) = I and 
Di (0) = Pi'. Hence the (unnormalised) state 

N 

|V'(d)):-(g)A(5)|AKLT), (7) 

i=l 

interpolates smoothly between |AKLT) at d = 1 and |G) 
at (5 = 0. If we restrict to the domain d > 0, then Di{6) is 
invertible. Given that |AKLT) is the unique state anni¬ 
hilated by every ^ follows that |?/'(<5)) is the 

unique state annihilated by every deformed interaction 
term 

[A(<5)-' G D,{5)-^] ® D,{6)-^] . 

( 8 ) 

We remark that these couplings take the form ABA 
where A and B are positive operators, thus the cou¬ 
plings themselves are positive. (It is not a simple unitary 
change of basis, which would leave the spectrum of the 
Hamiltonian unchanged.) Note that these couplings have 
an unbounded spectrum in the region d > 0 with some 
eigenvalues tending to infinity as d —> 0. 

To ensure that the spectral properties of the family of 
models is not obscured by the diverging spectral behavior 
of the local couplings, we replace the interaction terms 
with a set of projectors {Q{6)ij : {i,j) & E} where Q{6)ij 
is defined as the projection onto the image of Q. 

Thus we define a family of two-body Hamiltonians 

H{6) = ^ Q{S)., , (9) 

{id) 

which has \ip{6)) as a unique ground state. Given that 
the operators Q{S)ij are well-defined for any d > 0, and 
that lim 5 _>o \tp{6)) = |G), we conclude that the ground 
state of H{S) can be made arbitrarily close (in terms of 
fidelity) to |G) by setting d > 0 sufficiently small. 

The interaction terms Q{6)ij are defined for any d > 0, 
however at <5 = 0 the operators D[5) become singular 
and thus Eq. (|^ is not defined. We may neverthe¬ 
less define the Hamiltonian at <5 = 0 by taking the limit 
(5(0) := lim5^oQ(i^) and letting iL(0) := Q(0)q- 

This limit exists, and can in principle be computed by 
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applying Gram-Schmidt orthogonalisation to the column 
span of Eq (|^ for non-zero S, then taking the limit as 
d —>■ 0. H{S) is continuous and Q{S) have the same rank 
in the entire interval S € [0,1], however we note that the 
ground state degeneracy of ff(S) increases drastically at 
3 = 0. This behaviour is studied in greater detail in 
the following section. This essentially implies that the 
ground state of ff(3) cannot be made exactly equal to a 
graph state by setting (5 = 0, despite the fact that it can 
be made arbitrarily close by setting 3 small and non-zero. 


III. SPECTRAL PROPERTIES 

The construction detailed in the previous section pro¬ 
vides us with a family of two-body, frustration-free anti¬ 
ferromagnetic Hamiltonians ff(3) of Eq. Q that interpo¬ 
lates between a Hamiltonian possessing an AKLT state 
as its unique ground state for i5 = 1 and a Hamiltonian 
possessing a unique ground state that is arbitrarily close 
to a graph state for d —)■ 0. Like the AKLT Hamilto¬ 
nian neighbouring interaction terms Q{3)ij do not 
commute, and analysing spectral properties of H{3) is 
non-trivial. In this section, we study spectral properties 
of the Hamiltonian H{3) of Eq. by making use of 
the spectral properties of the simpler Hamiltonian H'^ of 
Eq. (|3|. 

In particular, we will study the behaviour of the spec¬ 
tral gap separating the ground state energy and the first 
excited state. To be precise, for hxed integer N, we will 
define the gap Ajv = Ei — Eq to he the difference be¬ 
tween the smallest and second smallest eigenvalues of the 
Hamiltonian with system size N. Note that Eq = 0 due 
to our Hamiltonian being a frustration-free sum of projec¬ 
tors. We say that a Hamiltonian is gapped if there exists 
a constant A > 0 and integer M such that A^r > A for 
all N > M. 

Proving that a model is gapped is hard in general [SD] . 
Whether or not the AKLT Hamiltonain E[^ is gapped on 
two-dimensional lattices remains an open question. The 
family E[{3) of Hamiltonians contains the AKLT Hamil¬ 
tonian iL(1) = and thus proving that H{3) is gapped 
for all (5 > 0 is at least as hard as proving a gap for E[^. 
We use a range of methods to study the spectral gap in 
the region where 3 is small and the ground state is close 
to a graph state. 

Note that there are two key parameters determining 
the size of the gap of H{3): the parameter 3 (which re¬ 
flects the closeness of the ground state to a graph state), 
and the system size N. First we will consider the case 
where N is fixed and 3 is allowed to vary. We will show 
that the gap, as a function of 3, approaches zero as 3 ap¬ 
proaches zero. In other words, there is a tradeoff between 
the fidelity of the ground state with a graph state, and 
the gap of the Hamiltonian. 

While small 3 implies small gap for finite-sized sys¬ 
tems, we will then show in Sec. Ell that small 3 also 
guarantees the existence of a gap in the thermodynamic 


(inhnite system size) limit. In other words, we prove that 
for 3 hxed and sufficiently small, there exists a constant 
A independent of system size, for which Ajv > A for all 

N. 


A. Behaviour of the gap as <5 —> 0 for finite lattices 

Here we will prove that, for hxed system size, there is 
a tradeoff between the gap of the Hamiltonian and the 
hdelity of the ground state with a graph state. We stated 
this result previously in [53], although did not provide a 
rigorous proof. 

The Hamiltonians we consider are not related simply 
to the stabilizer Hamiltonian for a graph state (which is 
clearly gapped). In fact, there is no point in the family 
where the ground state is exactly a graph state. We 
show in this section that for any hnite sized system the 
gap shrinks to zero in the limit where the fidelity of the 
ground state with a graph state tends to one. This result 
necessarily comes about because we have restricted to 
two-body Hamiltonians for which certain ground states 
(e.g., graph states) cannot be reached without closing the 
gap. 

We will first prove a result about a class of graph states 
which we call non-trivial graph states. A graph state is 
non-trivial if the graph G on which it is defined satisfies 
the following two properties. Firstly, the degree of each 
vertex of G is at least two. Secondly, no two distinct 
vertices in G have an identical set of neighbours. These 
properties are satisfied on all 2D lattices considered in 
this paper. We have the following lemma: 

Lemma 1. Let |G) he a non-trivial graph state on a 
graph of spin-112 particles. Let A := {i,j} he any two 
sites and B he its complement. Then the reduced density 
operator pA — trB(|G)(G|) is \l. Moreover, the only 
positive, two-hody, spin-1 f2 Hamiltonian with |G) as a 
zero energy eigenstate is the zero Hamiltonian. 

Proof. We will use the fact that the reduced density op¬ 
erator Pa of a stabilizer state with stabilizer S is given 
by 

(tGSa 

where Sa Q S is the set of stabilizer elements a G S 
which act as the identity outside A. Eq. ( |Io| ) will hold 
for any set of particles A, however here we will set A to 
be an arbitrary two-particle region {i,j}. 

Any element in the stabilizer a G S can be expressed 
uniquely as cr = n;! K'f '-, where Ki are the graph state 
stabilizer operators (defined in Eq. 0 and 7 / S {0,1} 
for each 1. It is clear that any a G S with 7 ; = I for 
I ^ A cannot be an element of Sa- This eliminates all 
but four possible stabilizer elements from membership in 
Sa- Ki, Kj, KiKj and I. As we have assumed that 
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the degree of G is at least 2, Ki and Kj must act non- 
trivially outside A, and thus Ki and Kj are not contained 
in Sa- Furthermore, from the assumption that i and j 
do not have an identical set of neighbours outside A, it 
is clear that KiKj is also not contained in Sa- There¬ 
fore Sa = {-f} and from Eq. (10), pA = \l- Thus any 
positive non-zero interaction term acting on A will have 
positive, non-zero energy. Since this holds for every pair 
of particles in the state, the only two-body Hamiltonian 
with positive interaction terms and zero energy will be 
the zero Hamiltonian. □ 


The above result holds for spin-1/2 Hamiltonians, how¬ 
ever it does not directly apply to our family of Hamilto¬ 
nians, which generally involve higher-dimensional parti¬ 
cles. We will now generalise to the case where each graph 
state qubit is encoded in a two-dimensional subspace, 
called the logical space "H/, of a higher-dimensional phys¬ 
ical particle. We find that any two-body, frustration-free 
Hamiltonian with such a graph state in its ground space 
must act trivially on the logical space of each particle. 

Lemma 2. Consider a system of N particles 
where the dimension di = 2j + \ of each spin-j system 
'W is at least 2. Let |C?) be a non-trivial graph state with 
N vertices defined such that the i-th qubit is encoded into 
a two-dimensional subspace of the i-th physical particle 
Hi ^ Let hij be Oj positive, two-body opevotov act¬ 
ing on W 0 W and satisfying hij\G) = 0 . Then hij is 
supported entirely on {Hf 0 Hf)^, i.e. hij = P^hijP^ 
where P^j is the projector onto {Hf 0 Hf)^ . 

Proof. Let A = {i,j} be a two particle region. Let Pl be 
the projection onto Hf 0 H'/. Using a decomposition of 
the identity I = Pl -\- P/-, we have 

hij = PLhijPL + {P^hijPL + h.c.) -H P^KjPl ■ (H) 

From hij\G) = 0 we have that PLhij\G) = 0 and 
Pijhij\G) = 0 . From Lemma we also have that the 
reduced density operator on the region A satisfies pA = 
\Pl- By rearranging tv{PLhij\G){G\hijPL) = 0 and 
tr(P/-hij|G)(G|%P/-) = 0, we obtain ||Pl%Pi,||2 = 0 
and ||P/'hijPi ||2 = 0 where ||P ||2 := ■\/tr(P, B) is the 
Frobenius norm. Hence only the last term in Eq. 0 is 
non-zero. □ 


This implies that any frustration-free, two body Hamil¬ 
tonian acting on qudits with a graph state as an ex¬ 
act ground state has the undesirable property that the 
ground space contains the entire logical space 
Thus the ground space of the Hamiltonian is (at least) 
2^-fold degenerate. From this it follows trivially that 
our Hamiltonian H[6) defined in Eq. © acting on a 
system of fixed size N with a non-trival graph state |G) 
as an approximate ground state, has a gap that shrinks 
to zero as 5 tends to zero. We prove the result generi- 
cally for any two-body, frustration-free Hamiltonian with 
a approximate graph state as a unique ground state. 


Theorem 1. Let H{6) = J2{i j) ® two-body, 

frustration-free Hamiltonian acting on N qudits, defined 
on the interval S € [0,1], where each term hij (6) is a 
projector that varies continuously in 5. Assume that for 
(5 > 0, H{d) has a unique ground state |V'(<5)) which has 
the property that lim^-^o If’id)) = |G), where |G) is a 
non-trivial graph state with N vertices (in the sense of 
Lemma^. Then lim^-^o Ajv = 0. 

Proof. By continuiW, each h{0)ij has |G) as a zero eigen¬ 
state. By Lemma ^ H{0) contains the 2^-dimensional 
logical space il® ground space. The fact 

that H{5) has a unique ground state while P(0) has a 
groundspace that is at least 2^-dimensional, and the fact 
that H{5) —>■ H{0) as 6 -i- 0, implies a large number of 
zero-energy crossings at 6 = 0. The existence of ground 
state energy crossings implies lim 5 _>o A^r = 0. □ 


We have thus shown for any finite-sized system, that 
if H (6) is a two-body, frustration-free Hamiltonian with 
a ground state that tends to a graph state as S tends to 
zero, then the gap of the Hamiltonian must also tend to 
zero as S tends to zero. We remark that if the bodiness 
of the Hamiltonian is greater than two, this trade-off is 
not necessarily observed. For instance, it was shown in 
Ref. [24] that for three-body parent Hamiltonians in one 
dimension, the fidelity can be improved arbitrarily with¬ 
out decreasing the size of the gap. 


B. Gap in the thermodynamic limit 

We now prove that the Hamiltonian H(S) is gapped 
in the thermodynamic limit if J > 0 is independent of 
N and sufficiently small. This statement may seem at 
odds with the result of the previous section, where we 
showed that the gap of the finite-sized system tends to 
zero as S goes to zero. The two results can be reconciled 
by the fact that we are taking two different limits: here 
we are considering the case where 6 is fixed and IV tends 
to infinity, while in the previous section we considered 
the case where N was fixed and S tends to zero. 

Our proof relies on properties of parent Hamiltonians 
for injective tensor network states. Injectivity is a generic 
property of tensor network states on a lattice ED, and is 
useful in our proof as it implies the existence of a unique 
ground state. To prove injectivity, we first perform a 
blocking and construct a blocked parent Hamiltonian, 
for which we can prove a spectral gap that is stable for 
small perturbations in 6 around zero. After defining ten¬ 
sor network states and blocked parent Hamiltonians, we 
then prove that the class of two-body antiferromagnetic 
Hamiltonians H{S) is gapped on various lattices for suf¬ 
ficiently small <5 and provide lower bounds on the size of 
the gapped region for trivalent lattices (i.e., the largest S 
for which we can prove H{6) is gapped). 
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1. Tensor network states and PEPS 

The projected entangled-pair state (PEPS) framework 
provides an efficient description of a class of multipartite 
quantum states in arbitrary spatial dimensions. Here we 
provide a brief description of the PEPS framework, as it 
applies to our problem; for further details, see Ref. [55]. 

A tensor is a multi-index array of complex numbers 
Aa^^a 2 ,...,an where each index has some finite dimen¬ 
sion dfc, which we call the bond dimension. Graphically, 
we represent a tensor as in Fig. Ba). In this graphi¬ 
cal framework, when two tensors are connected by an 
edge, their corresponding indices are contracted. For the 
purposes of defining tensor network states, we identify 
two distinct types of indices. The first type are physical 
indices, denoted /3i, associated with an orthonormal ba¬ 
sis {\Pi),i = l,...d} for each elementary d-dimensional 
quantum system. The second type are virtual indices, 
denoted using a indices, which are not associated with 
physical degrees of freedom but are contracted internally 
within the tensor network to allow for a non-product en¬ 
tanglement structure. For clarity, in our tensor network 
diagrams we will not usually draw an edge for each phys¬ 
ical index, instead implying the existence of a physical 
index on a tensor with a superscript /3 on the tensor la¬ 
bel. 

Consider a quantum state defined in the above ba¬ 
sis as \lp) = A 

projected entangled pair state (PEPS) is a state where 
is obtained by a contraction of a network of 
tensors, one tensor for each physical index, as in Fig. Sa). 
Graph states and AKLT states are well studied examples 
of PEPS. Appendix [A] presents a specific standard choice 
of tensors for describing graph states and AKLT states. 

2. Blocking the Hamiltonian 

When studying Hamiltonians that have PEPS as 
ground states, it is often simpler to consider a Hamil¬ 
tonian which has the same ground state as the physical 
Hamiltonian, but with interaction terms that act on a 
larger number of particles. Such Hamiltonians, which 
we will call blocked Hamiltonians, share many properties 
with the physical Hamiltonian, but are often simpler to 
study because the spectrum of the local terms can be cho¬ 
sen to be very simple. Here, we define a blocked version 
H^{6) of the parent Hamiltonian H{6). This blocking 
will assist in our proof of a gap for H{6), as we will show 
that proving a gap for H(d) is equivalent to proving a 
gap for H^{S) for (5 > 0. The blocked Hamiltonian will 
possess an injectivity property that allows us to prove a 
gap. In addition, H^{S) can be rigorously proven to be 
gapped not only for (5 = 0 but for small <5 > 0. 

We define the blocked Hamiltonian H^{S) as follows. 
Let G be a graph, and consider the PEPS describing the 
state \ip{S)) of Eq. 0 on this graph that interpolates 
between the graph state |G) at 6 = 0 and the AKLT 



(b) 


FIG. I: A connected region of tensors can be ‘blocked’ 
to form a single tensor by contracting all internal 
indices within the region, (a) Tensors in a five-particle 
region, (b) Blocked tensor. 

state at S = 1. Given a connected region R C V of n 
vertices and r outgoing edges we define the block tensor 
as the tensor obtained by contracting 
all virtual indices of tensors within the region and leaving 
outgoing virtual indices uncontracted, as illustrated in 
Fig. Qb). This blocking naturally defines a map from 
virtual, to physical degrees of freedom 

H(«)(<5)= ^ Hf(;;::;£:(R,5)|/3i,...,/3„)(ai,...,a.|. 

CKi. ,0:^ 

( 12 ) 

We will say that a region R of particles is injective if the 
associated map (S) is injective. Generally, if the ver¬ 
tices of a tensor network state |'0) can be partitioned into 
M disjoint connected regions V = U^^Ra such that each 
set Ra is injective, we say that the state |^) is injective. 

In Appendix we show that a region R for the state 
\tpiS)) is injective if and only if each vertex i € R has 
at most one outgoing edge (i.e., an edge connected to a 
vertex not in R). We will say that {Ra : a = 1,... ,M} 
is an injective covering of V if Ra is injective for every a, 
RaPRb — 0 for any a ^ b, and U^^Ra = V. Examples of 
lattices with injective coverings are illustrated in Fig. 
with the injective region highlighted. 

Let Ea^b Q E he the set of edges in G that connect a 
vertex in Ra to a vertex in Rb- Given a partitioning of 
the graph into injective regions {Ra : a = I,..., M}, we 
define a coarse-grained graph G' = {V',E') as follows. 
Each vertex in V' corresponds to a region Ra, and two 
vertices in V, corresponding to the regions Ra and Rb, 
are connected by an edge in E' if and only if Eafi is non¬ 
empty. 

Given the state {tpiS)) defined on a graph with a injec¬ 
tive covering, let .Bfo’^^((5) be the injective block tensor 
for the region RaURb, where Ra and Rb are neighbouring 
injective regions (where ‘neighbouring’ means that they 
are connected in G'). 

We define the blocked parent Hamiltonian for \'tp{5)) 
as 

H^{5):= nfo’'')(<5), 

{a,b)eE' 


( 13 ) 













7 




FIG. 2: The star lattice (a) and honeycomb lattice (b) 
with injective regions circled. A region is injective if 
and only if each particle has at most one outgoing edge 
(i.e. connected to a particle outside the region). 
Furthermore, we say these lattices are coverable as both 
can be decomposed as as a disjoint union of such 
regions. 


arbitrary graphs) we will impose some further simplifica¬ 
tions. We require that there exists an injective covering 
for the graph {Ra : a = 1,..., M}, for which both the 
number outgoing edges r and the number of particles n 
in each Ra is a constant independent of a and the system 
size N and that the block tensor for each Ra are 

all identical up to a unitary acting on the output for all 
(5 > 0. It would vastly complicate the proof if the injec¬ 
tive covering varied with the system size or if the spectral 
properties of the block tensor were different for different 
regions. 

We will call a graph satisfying all of the above proper¬ 
ties coverable. All trivalent Archimedian lattices are cov¬ 
erable: the honeycomb (6^), the square octagon (4,8^), 
the cross (4,6,12) and the star (3,12^). However, the 
square lattice (which does not have a injective covering), 
is not coverable. We will discuss generalisations non- 
coverable lattices in IIIIH51 


where n^“’^)((5) is the projector onto the orthogonal 
complement of the image of and the sum is 

taken over all neighbouring pairs of regions Ra and Ri,. 
Note that this Hamiltonian is not two-body: each term 
n(“’^)((5) will in general act non-trivially on all particles 
in the region Raii Rt- Given that the reduced^ density 
operator of [■(/'(iJ)) is supported on the image of B^‘^’^\S), 
it is clear that \ip{S)) is a ground state of H^{6). As is 
shown in Ref. |31| the injectivity condition further im¬ 
plies that |■^/’(<5)) is the unique ground state of H^{5). 

We thus have two Hamiltonians, H{5) and H^{S), 
both possessing \tp{6)) as a unique ground state. The 
Hamiltonian H{S) is two-body while H^{S) is not. A 
crucial difference between these two Hamiltonians is that 
H^{6) is well defined at 5 = 0 with a unique ground state 
exactly equal to the target graph state |G). On the other 
hand, as we showed in Sec. |HI A[ the two-body Hamilto¬ 
nian H{S) has an exponentially degenerate ground space 
at 5 = 0. In the following section, we will use the fact 
that H^{0) has a unique ground state and is gapped to 
prove that H{S) is gapped for sufficiently small 6. 


3. Restrictions on interaction graph 

We have imposed a number of constraints on the inter¬ 
action graph of the Hamiltonian, which we briefly sum¬ 
marise here. First, we have restricted to three-colourable 
graphs, for which the methods of Sec. |H B| can be used 
to project the AKLT state locally to a graph state on 
the same graph. To prove a gap in the thermodynamic 
limit, we will further require injectivity for our blocked 
Hamiltonian. As described in the previous section we can 
ensure injectivity by imposing that our graphs have an 
injective covering. 

As we are primarily interested in lattices (rather than 


4. Proof of gap 


Here, we prove the following theorem regarding the 
existence of a spectral gap for H{6), which is the central 
result of this section. 


Theorem 2. Let H{d) be the two-body Hamiltonian de¬ 
fined in Eq. ([^ on a coverable graph. Then there exists 
(5c > 0 such that H{S) is gapped for all S in the open 
interval (0,5c)- 

The proof of Theorem [^follows from three simple lem¬ 
mas. We will first show that the blocked Hamiltonian is 
gapped at 5 = 0. 

Lemma 3. H^{Q) is gapped. 


Proof. We prove this lemma using the fact that |'0(O)) = 
|G), and by showing that this graph state can be locally 
transformed to a product of Bell pairs via a product of 
unitaries that act only within injective regions (and not 
between them). Under this unitary map, the Hamilto¬ 
nian can be seen to be trivially diagonalisable, commut¬ 
ing and thus gapped. 

Let G be a coverable graph an d co nsider its tensor 
description, provided in Appendix A 2 Let {Ra : a = 
1,...,M} be an injective covering of G. Let Wa ■= 
0^ - 0 T)ij be a product of CZ gates acting 

on all pairs of neighboring particles in the region Ra 
where each {CZ © T)ij acts as CZ on the logical sub¬ 
space Hf 0 Rj C "Hi 0 "H of particles i and j and as 
the identity on the remainder of the space. Using the 
identity given in Eq. (A21 on the tensor description of 
H(“’^)(5) on two injective regions Ra and Rt, we see that 
the image of [Wa 0 Wb]B^°'’^\0) becomes a simple tensor 
product. To be precise, let denote the interior of the 
region R, i.e., the set of vertices in R that are not con¬ 
nected to vertices outside R, and let OR := R\R^ denote 










(15) 


the boundary of R. The image of [Wa 0 (0) is 


1 +)* 


\H)jk 


n 




. i, 


a&d(RaURb) 


(14) 

where \H)jk G Rf 0 is a two-qubit graph state on 
particles j and k and Ea^b was as previously defined as 
the set of edges in G that connect region to region 
Rb- 

Let be the projector onto the space orthogo¬ 

nal to the space defined by Eq. (14). The product form 


of Eq. ( |14[ ) immediately implies that the set of opera¬ 
tors : (a, b) G E'} pairwise commute. It is also 

clear that the Hamiltonian H' := 'J2{a b)eE' has a 

trivial unique ground state consisting of \E[) on all con¬ 
nected pairs of vertices that are from separate regions, 
and |-|-) on the remaining vertices that are only connected 
to other vertices within the same region. Hence El' is 
gapped. Einally, given that H' and H^{0) are related by 
conjugation by Wa we conclude that H^{0) is also 

gapped. □ 


Having shown that H^{0) is gapped, we next show 
that H^{5) is gapped for a finite interval (0,5c)- This 
result can be seen to follow quite directly from stability 
results for canonical parent Hamiltonians [32], but we 
provide a self contained proof of this fact in Appendix [Pj 
based on methods used previously in Refs. [331 [33] for 
frustration-free Hamiltonians. The idea of the proof is to 
show that a PEPS parent Hamiltonian is gapped when 
its ground state is sufficiently close to an isometric PEPS 
(as is the case when 5 is sufficiently small). 

Lemma 4. There exists 6c > 0 such that H^{6) is 
gapped for all 6 in the interval [0, 5c). 

Proof. See Appendix]^ □ 


Thus, there is an open interval on which we can prove 
that the blocked Hamiltonian H^{6) is gapped. To com¬ 
plete the proof of Theorem]^ we need only to show that 
if H^{S) is gapped for some 5 > 0, then the two-body 
Hamiltonian H{6) is also gapped. Here, and for the rest 
of the paper, we will use the convention that if A and B 
are Hermitian operators, then A> B means that A — B 
is a positive operator (i.e., each eigenvalue is real and 
greater than or equal to zero). 

Lemma 5. For 5 > 0, H{6) is gapped if and only if 
H^{6) is gapped. 

Proof. Let {Ra : a = 1,..., M} be an injective covering 
of G such that the degree of each coarse-grained vertex 
in V (i.e., the number of outgoing edges of each Ra) 
is a constant r (the existence of such a covering was an 
assumption of our graph). We rewrite H(5) by grouping 
interaction terms according to the coarse-grained graph. 


as 


H{6)= Y. 

{a,b)GE' 

where the grouped terms can be expressed in terms of 
this new course graining as 

h{S)ab-.= Y Y > ( 16 ) 

{i,j)eEaUEi, {i,j)eEab 


where C if is the set of edges contained entirely in the 
region Ra and, as before, if^.b is the set of edges joining 
vertices in Ra and Rb. It is straightforward to show that 
the sum in Eq. (15) is equal to H (5) as defined in Eq. (|^ . 

Comparing h{S)ab to the blocked interaction term 
n(“’'')(5), one finds that the kernel of a single term h{S)ab 
is equal to the kernel of the projector n(“’^)(5). We will 
not provide a detailed proof here; it is possible to ver¬ 
ify following the inductive proof in Ref. [25j . Appendix 
A, where it was shown that the spin-3/2 AKLT model 
(and deformed versions of it) has a unique ground state. 
Because these operators have the same kernel, we have 

A„,i„(5)n(“’^)(5) < h{S)ab < A„,ax(5)n(“’^)(5), (17) 


where Ainin(5) > 0 and Aniax(5) > 0 are respectively the 
smallest and largest non-zero eigenvalue of h{6)ab- Im¬ 
portantly, Amin (5) > 0 and Amax(5) > 0 are indepen¬ 
dent of the system size. Erom Eq. 0 , it follows that 
Amin(5)i7^(5) < H{6). Therefore, for any system size N, 
if H^{6) has gap A(5), H{5) will have a gap of at least 
Amin(5)A(5). Hence for any fixed 5 > 0, H{5) is gapped 
if H^{S) is gapped. Likewise, as X^a^{d)H{6) < H{S)^ , 
E[{6) being gapped implies H^{5) is gapped. 

We remark that, while Amin(5) > 0 for all 5 > 0, taking 
the limit gives lim^-^o Amin(5) = 0. Thus we cannot prove 
a gap at 5 = 0 (which is unsurprising, given that we 
showed that the gro und space becomes degenerate at this 
point in Sec. HI A[ ). □ 


This completes the proof of Theorem that, under 
certain assumptions on the interaction graph, H{5) is 
gapped for 5 sufficiently small. 


5 . Generalisation to arbitrary lattices 


In the proof of Theorem we showed that the two 
body Hamiltonian H (5) is gapped for sufficiently small 
5 > 0 provided the int eraction graph is coverable (as 
described in Sec. |HIB~3 ). Here we will discuss obstacles 
to generalising this result to non-coverable graphs. 

Consider, for example, the square lattice. The state 
|'!/'(5)) is the unique ground state of the two-body, spin- 
2 Hamiltonian H{S) for 5 > 0 on this lattice, but is it 
gapped? The state |V'(5)) is not injective on this lattice, 
and thus Theorem]^ does not apply to it directly. We can 
nevertheless follow the proof of Theorem]^ and construct 
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a parent Hamiltonian for every 5 > 0 with \iIj{S)) as a 
unique ground state. 

First we partition the particles into 2x2 squares {Ra : 
0 = 1,... ,M}. Let be a blocked interaction 

term that acts on a 2 x 2 square of four such regions Ra, 
Rb, Rc, and Rd, which is defined, as usual, as the projec¬ 
tion onto the orthogonal complement of the image of the 
block tensor show that the 

blocked Hamiltonian H^{d) := J2{a b c 
with sum taken over all such squares, has \iIj{S)) as a 
unique ground state for both (5 = 0 and (5 > 0. It is 
also possible to show that Lemma|^holds (i.e., H^{0) is 
gapped) and Lemmaholds (i.e., the two-body Hamil¬ 
tonian H{d) is gapped if and only if H^{S) is gapped 
for (5 > 0). The problem arises when we try to interpret 
H^{6) for (5 > 0 as a small perturbation to the gapped 
Hamiltonian iJ^(O), as is required to prove Lemma ^ 
One finds that there is a discontinuity of the rank of the 
interaction terms n(“’^’'^’‘^)((5) at d = 0, specifically the 
rank is smaller at (5 > 0 than at (5 = 0. This is in con¬ 
trast to the injective case, where the rank of the block 
tensor (and therefore also the interaction terms) is con¬ 
stant for all S > 0. As a result, the stability arguments 
used to prove Lemma and the proof methods of Ap¬ 
pendix cannot be applied directly. Hence, it does not 
seem as straight-forward to prove a gap for non-injective 
lattices, as for injective ones. 


C. Gapped regions for trivalent lattices 


Theorem states that for the two-body Hamiltonian 
ff(S) (defined on an appropriate lattice) there exists a 
(5c > 0 such that R(S) is gapped for all (5 on the open 
interval (0, dc)- Here we will explicitly compute lower 
bounds on Sc on various trivalent lattices. 

Let R(S) be defined on a graph with an injective cover¬ 
ing {Ra : a = 1,...,M}, where each injective region Ra 
has the same number of outgoing edges r. Consider the 
operator {S)B^^‘‘'> {S) and let 7min(<5) and 7max(<5) 

be its smallest and largest eigenvalues respectively. Note 
that the ratio 7max(<5)/7min(^) is a measure of how far 
jg f]-om the identity. At d = 0, this 
ratio achieves its minimum of 1 (i.e., i?(^“ii(0)i?(^“i(0) 
is the identity). In Appendix we show that H{d) is 
gapped if 

Tmax (<5)/7min ((5) < /xq, where 


I I 
- 2 2 



(18) 


Thus our bound depends on the number of outgoing 
edges in each injective region r. 

Note that /xq > I is greater than I for any r, and 
thus there is always a region around (5 = 0 for which this 
bound proves that the Hamiltonian is gapped. 

We can evaluate the operator {S)B^^‘^^S) 

by contracting a tensor network, illustrated in 
Fig. For r = 3,4,5,6 we have « 



A 

iY - 

A 

iY - 

A - 

iY 


(Sy D 

2(Sf D 

3(<5) 

2 


A 

iY - 

A 

iY - 

A - 

.y 




FIG. 3: Expressing the operator B^^<‘'> {S)^{S) for 
the star lattice (with r = 3) as a contraction of tensors. 
Where the indices at the top represent the input and 
the indices at the bottom represent the output. 
Physical indices in this figure are represented by red 
edges. The Hamiltonian H{S) is gapped when the 
eigenvalues of this operator satisfy Eq. (181. 


0.828, 0.873, 0.899, 0.916 respectively corresponding re¬ 
spectively to Sc « 0.28, 0.13, 0.12, 0.08 and the Hamil¬ 
tonian H{S) will be gapped for all S G (0,(5c). Our upper 
bound on the size of the gapped region decreases as a 
function of r, the number of outgoing edges in an injec¬ 
tive region. 


IV. CONCLUSION 

We have shown that methods from quantum informa¬ 
tion theory, in particular from the study of graph states 
and AKLT states for quantum computation, can be use¬ 
ful for proving spectral properties of 2D antiferromag¬ 
netic spin models. In particular, for a family of 2D anti¬ 
ferromagnetic spin models obtained by a one-parameter 
deformation of the 2D AKLT model defined here, we 
have proven that there exists a spectral gap separating 
the ground state from the first excited state in the ther¬ 
modynamic limit for a range of parameter space of this 
family. 

An implication of our results is that this family of mod¬ 
els has many of the desired properties for quantum com¬ 
putation, in particular that there exists a finite region in 
the parameter space that is both computationally uni¬ 
versal and gapped. The existence of a gap may allow 
for efficient methods to prepare such resource states for 
quantum computation [55H57] . 

It would be desirable to prove that the gapped re¬ 
gion extends to (5 = 1 for the class of lattices defined 
here, where ff(S=l) is the AKLT Hamiltonian This 
would settle the long-standing conjecture regarding the 
existence of a spectral gap for the spin-3/2 2D AKLT 
model. While the region we prove to be gapped unfortu¬ 
nately does not extend to (5 = 1, there are various ways 
by which proofs of a larger gapped region may be sought. 
For instance in, mixing properties of Markov chains were 
used in Ref. [31] to show that the Ising PEPS parent 
Hamiltonian is gapped up to the critical point. Other 
methods from statistical physics, e.g., using ‘martingales’ 
[55] may also be used to improve our result. 
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Appendix A: PEPS definitions 


In this section we will provide simple tensor-network 
descriptions of AKLT states and graph states and for the 
interpolating path between the two. 


1. AKLT states as PEPS 


An AKLT state may be defined on a graph G = {V, E) 
as follows. At each vertex i of the graph place a spin-S'^ 
particle, where Si = —1)/2 and di is the degree of ver¬ 

tex i G V. A basis for the Hilbert space of each particle 
is given by {\M) : M = Si, Si — 1,, —Si + 1, —Si}, cor¬ 
responding to the eigenstates of the z-component of the 
spin operator S^ . For each vertex i let be a 

tensor of Clebsh-Gordan coefficients for a set of di spin- 
1/2 particles coupling to total spin-S'i, where fdi = M is 
the total spin along the z-axis and ak is the component 
of the fc-th constituent spin-1/2 along the z-axis. For 
example, if Si = 3/2, aa.aa have the following 
non-zero entries 


Graph states may also be represented as simple tensor 
networks. Let G = {V,E)he a. graph defining a graph 
state, according to Eq. Q and let di be the degree of 
vertex i. At each vertex i G V place a spin-1/2 particle 
and define a tensor G with di virtual indices and for which 
there are only two non-zero entries 

^ 00...0 = ^ 11...1 — 1 ) ( A . 1 ) 


where the labels 0 and 1 for physical and virtual indices 
may be regarded as labels for spin 1/2 and —1/2 states 
respectively. Applying a GZ gate to the physical index 
of any two disconnected G tensors is equivalent to adding 
an additional virtual index to each G and contracting a 
Hadamard matrix El = -^(X + Z) between them. 

Graphically this can be represented as 



where CZ^J’o? 

Pi >P2 


is the tensor of the CZ in the standard 


computational basis with /3j, being input indices and 
Pi, 132 being output. Although we have illustrated this 
with with C tensors that initially have three virtual in¬ 
dices, this identity will hold for any number of initial 
virtual indices (including zero). Note that as = H, 
there is no need to specify which index of H is contracted 
with which neighbour. The tensor network for a general 
graph state therefore consists of a C tensor at each ver¬ 
tex and a H contracted between neighbouring G tensors 
on each edge, as illustrated in Fig. |^b). 


3. The family of states IPiS)) as PEPS 


^000 = — 1 ! 

Aqoi = Aqj^q = A/qq = 1/-y/S, 

Aqii = ^110 ~ ^101 ~ 1 /1 

where 0 and I subscripts refer to the spin eigenstates 
rUz = 1/2 and niz = —1/2 respectively. 

Given the A tensors, the resulting tensor network for 
the AKLT state is defined according to Fig. |^a). An lY 
matrix is placed on each edge connecting neighbouring A 
tensors. This is due to the fact that the bonds used to 
define the AKLT PEPS are singlets l/-\/2(l C) fF)(|00) -b 
|11)), rather than l/-\/2(|00)-I-111)), and gives the AKLT 
state its antiferromagnetic character. (Note that we do 
not need to specify which index of these edge lY tensors is 
contracted with which neighbouring tensor, as swapping 
these indices only results in an overall ignorable phase 
of —1.) Note that we may contract tensors without a 
physical index (the lY tensors in this case) with their 
neighbours such that each tensor has one physical index, 
thus satisfying the usual definition of a PEPS. 


Finally, we will define the tensor network for the family 
of state |V'(<^)) defined in Sec. II B| interpolating between 
between the AKLT state (at 6 = 1) and a graph state 
(at ^ = 0). As these states can be obtained by apply¬ 
ing a product of single-particle operators ®iDi{5) to 
the AKLT state, the tensor network of \ip{6)) can be ob¬ 
tained simply by replacing the A tensors of the AKLT 
state with the (5-dependent tensors A((5)^/_^2,...ad. ■= 

. We show in Appendix B that 
the tensor network for the state I'PiO)) is equivalent, up 
to local unitaries, to the tensor network of a graph state. 

We remark that the tensors we have used to define 
AKLT states, graph states and |^/'((5)) have only real en¬ 
tries. 


Appendix B: Local conversion of AKLT states to 
graph states 

Here we will provide a proof of the fact that, on three- 
colourable graphs, AKLT states can be converted to 













11 



FIG. 4: PEPS for an AKLT state (a) and a graph state 
(b) defined on a particular graph of 6 vertices, where 
the A and C tensors are defined in the text. Note that 
the coloured tensors do not have a physical index, and 
can be absorbed into their neighbours if desired. 


graph states by applying a rank-2 projector to each par¬ 
ticle. The result was originally proved in Ref. [TO] for 
the particular case of spin-3/2 particles using the stabi¬ 
lizer formalism. Here we will provide an alternative proof 
using tensor networks. 


Let G be a three-colourable graph and let |AKLT) be 
the A KLT state defined on this graph, according to Ap¬ 
pendix A 1 Let {ci € {x,y, z} : i = 1... N} be a three- 
colouring of the vertices of G, in the sense that Ci ^ Cj if 
i and j are neighbours in G. 

Define and as in Eq. (^. Then the AKLT 

state can be transformed into a graph state by these local 
operators in the following sense. 



(a) (b) 



(c) (d) 


FIG. 5: Illustrating how to convert AKLT states to 
graph states with single-particle projectors. We 
consider an AKLT state on a triangle where a different 
projector is applied to each particle, (a) AKLT state on 
a triangle, (b) State after applying P“, P^, P^ to 
particles 1, 2 and 3 respectively, (c) Rewriting each 
edge matrix as H multiplied on the left and right by 
z-rotations. (d) Absorbing all z-rotations into the G’s. 
They will appear only on the physical indices (not 
pictured). 


Lemma 6. The following relation holds: 


are respectively 


N 


N 


|P"‘|AKLT) (x(g)Z(0O|G) 


(Bl) 


where |G) is the graph state as defined in Appendix A 2 


Z{9) = diag(I, exp(i0)) is a local z-rotation by 0 and 
where Oi G {0,7r/2, tt, 37r/4} V i. 


Proof. Applying P'^' to the AKLT state attaches 

one of P^, py, P^ to the physical index of each tensor. It 
is straightforward to show that applying P^ transforms 
the AKLT tensor A to the graph state tensor G, i.e. 

= Gf,Now let be 

the spatial rotation by angle 9 about the axis r, where 
S = {S^, , S^) is the vector of spin operators. (Note 
that we have adopted a different convention in our above 
definition of the z-rotation of a spin-1/2 particle Z{9), 
so that, for example Z{Tr) equals Pauli Z, however these 
two definitions differ only up to an overall phase.) From 
the symmetry of the A tensors, applying the same spa¬ 
tial rotation to the physical index and all virtual indices 
leaves A invariant. Furthermore, applying either P“ or 
P^ is equivalent to applying a rotation to the physical 
index, followed by applying P^. From this, one can see 
that, up to an overall phase, X)/ 3 ' 


E 

TT TT 







E 

H' H' 

^ Ct2-,Oi'2 








(B2) 


where H = -^{X + Z) and H' = + Z) are complex 

Hadamard matrices. Applying the projectors P”^* 
to I AKLT) thus converts all A tensors to G tensors and 
attaches H and H' to their virtual indices according to 
which projector (P^ P^ or P^) was applied. This is il¬ 
lustrated for a three particle AKLT state in Fig[^a)-(b). 
To prove the lemma, we show that this tensor network 
is equivalent to the tensor network of a graph state de¬ 
scribed in Appendix | A 2| up to local z-rotations. We will 
use the fact that for the G tensor, Z{9) applied to any 
virtual index has the same effect as Z{9) on the physical 
index, i.e., 


Z{9) 




= ^(^')/3,/3'Gf^„2 

(B3) 

for any k. 

As seen in Fig[^b), aft er a pplying the projectors Q to 
the state and using Eq. (B2), we are left with a product 
of matrices along edges of the tensor network between 
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the C tensors. By assumption, if i and j are neighbour¬ 
ing vertices, Ci ^ Cj. Thus there are only three distinct 
edges we must consider xy, xz and yz. These edges have 
the matrices iHYH', iHY, and iYH' applied on them 
respectively. These can all be expressed as H multiplied 
on the left and right by z-rotations and by an overall 
phase, specifically 


iHYH' = 

iHY = ZHZ , 

iYH' = iSHS , (B4) 


where S = Z(7r/2). Using (B3) all z-rotations can be 
moved to the physical indices of the C. From this we 
see that the tensor network is equivalent to the tensor 
network of a graph state described in Appendix |A 2[ up 
to z-rotations on the physical indices. □ 


We remark that from the definition of Di{5) in Eq. (§, 
we have P'''Di{S) oc and therefore the result also 
holds for all deformed states |'0(<5))- 


Corollary 1. Let G be a three-colourable graph, and 
{ci G {x,y,z} : i = 1... A^} be a three-colouring of G. 
Define the corresponding one-parameter family of states 
|'0(i5)) according to HB Then 


N 


N 


\p-'\fi{S))^^z{e,)\G), 


(B5) 


i=l 


i=l 


for all 5 > 0. 


Appendix C: Injective regions 


Here we will prove a necessary and sufficient condition 
for a region of particles in the tensor network of \fi{S)) 
to be injective. 

Lemma 7. A region R in the tensor network state 
for any S > 0 is injective if and only if each vertex in R 
is connected to at most one vertex outside R. 


Proof. Let i? be a region for which every vertex is con¬ 
nected to at most one vertex outside R, and b e 

block tensor for this region (as defined in Sec. IIIB2). 
First we will show show that this map is injective, i.e., 
R is injective. We showed in Appendix that we can 
transform the tensors of the state \tp{S)) to those of 
|^/>(0)) = |G) via a set of rank-2 projectors. Specifically, 


we have that (ZiieRPi) B^'^^S) oc B^^\0). Therefore the 
rank of B^^'>{6) is equal to or greater than the rank of 
5(fl)(0) for all i5 > 0. Thus showing B^^^O) is injective 
will imply injectivity of for all 5 > 0. 

To show that B^^^O) is injective, we will apply a uni¬ 
tary to its output (which will not change its injectivity) 
and transform it into a form that is clearly injective. We 
define a unitary W := 0^^ j)g£;(CZ 0 I)ij, as a product 



FIG. 6: Disentangling a region in the tensor network of 
a graph state by applying the unitary W. One can 
clearly see that the corresponding map from outgoing 
virtual indices to physical indices is injective, (a) A 
particular region in the tensor network of a graph state 
with three outgoing indices, (b) The same region after 
the unitary W has been applied to physical indices. 


of CZ gates acting on all pairs of neighboring particles 
in the region R where each (CZ 0 I)ij acts as CZ on 
the logical subspace Rf 0 'hj' C "Hi 0 "H of particles i 
and j and as the identity on the remainder of the space. 
(Given, however, that the image of B^^^ (0) is contained 
entirely in the logical space the action of W 

on the remainder of the space is not actually important.) 

We have illustrated the tensor network describiM the 
map B^^^O) before and after applying W in Fig. ^Us¬ 
ing Eq. ( |A2[ ) to remove edges, we see that W1^\0) 
describes the identity map from each outgoing virtual 
degree of freedom into the logical subspace Rf of the par¬ 
ticle i it is incident to, with all interior particles (those 
without an outgoing edge) placed in the |0) state. The 
identity map is clearly injective, therefore B^^^{d) is in¬ 
jective for any i5 > 0. 

It remains to show that having at most one outgoing 
index per particle is also a necessary condition for injec¬ 
tivity. Consider a region R containing a vertex with two 
outgoing indices ai and 02 , say. Let be the 

tensor corresponding to this vertex. From its definition, 
A is symmetric under any permutation of the virtual in¬ 
dices. Thus the input |01) — |10) on indices ai and a 2 
of the map will map to zero. Thus B^^\\) is 

not injective. As B^^\5) oc 0-^^^ Di((5)i?(^)(l) for all 
(5 > 0, (1) being non-injective implies B^^'> {S) is non- 

injective for all <5 > 0. □ 


Appendix D: Stability of gapped region 

Here we will provide a proof of Lemma which states 
that under the assumptions of Theorem]^ there exists 
(5c > 0 such that H^{S) is gapped for all 6 G [0,(5c). 

We make use of an argument, previously used in 
Refs. [SSI El] for frustration-free Hamiltonians. Let 
H = J2{jk) ^ two-body Hamiltonian acting on 

N particles, and for simplicity we will assume that the 
interaction graph has constant degree r. Assume that the 
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Hamiltonian is frustration-free and has a unique ground 
state with zero energy. If there exists a constant A > 0 
such that > AH for all system sizes, then H is 
gapped. 

Without loss of generality, we assume that the Hamil¬ 
tonian is scaled such that h\i > h^i for all interaction 
terms. A sufficient condition for the gap is that there 
exists a A > 0 such that for every particle c and for 
any pair of distinct interaction terms hie, hej that act 
non-trivially on c we have 

hichej ^ 2(r • (^1) 

To see this, we expand to get 


+ 

Vo'.fc) / 

(D2) 

The first term is greater than or equal to H. For the 
second term, we can split the sum into terms where hu 
and hpq overlap at one site, and terms where the inter¬ 
actions do not overlap (disjoint). The disjoint term is 
clearly positive, thus 


Ak,l),{p,q)-.{k,l)^{p,q) 


H^>H 


E hkihpq I . (D3) 

<{k^l)^{p^q) overlapping 


Using condition Eq. (Dl) in Eq. ( |D3[ ), we find that the 
right-hand side reduces to AH. Therefore H is gapped. 

Next, we show that if a Hamiltonian satisfying the 
gap property Eq. (Dl) is deformed by a sufficiently small 


amount, it will remain gapped. 

Lemma 8. Let H = ® frustration-free 

Hamiltonian with a unique ground state on a graph of 
constant degree r, suc h tha t are projectors satisfy¬ 

ing the gap condition (Dl) for some 0 < A < 1. Define 
the Hamiltonian H' = j) Kj with interaction terms 


— (A^ Aj) h^j (A^ Aj) 


(D4) 


for a set of invertible positive operators {Ai} which satisfy 
Ai > I and let p > 1 be the largest eigenvalue of A? 
(which we assume is independent of i). Define 

A'■.= !-p[l-A + 2{r-l){p-l)] . (D5) 

Then H' is gapped when 0 < A'. 

Proof. Note that A' < A for all p> 1. We claim that, 
1-A' 


^ij k T k 


' u/ 

ij 


2(r-l) 


{K,+hb^)>Q. (D6) 


Therefore if there exists p such that 0 < A', the deformed 
Hamiltonian will itself satisfy Eq. (Dll and be gapped. 
We define 0 = A^ — 1 > 0, which satisfies {p — 1)/ > 0. 


To prove Eq. (D6), multiply the left hand side on the left 
and right by A“^ 0 A“^ (g) A^^ and and simplify as follows 

1 — A' 

hijA^hjk + hjkA^hij -\- ® -^k ^ T ^ ® j 


- 2(r-l)' 

1-A' 


P hijAjhjji^ -)- hjj^Ajhij -)- 
— h2Ah.j'/i; -)- hjf^Ajhij 


2p{r - 1) 


{hij (g 7 -I- J (g hjk). 


Hj^'-j'T-jk 

1 


[1 - A + 2(r - l)(/r - 1)] {hij (g 7 -|- 7 g) hj/,) , 


2(r-l) 

— hijQjhjf^ -\- hjj^Qjhij 

hij hjk ~t“ hjj^hij -\- hjjf) 


-\- {p Xfihij -)- hji ^), 

P hij^jhjfz -(- hj].^Qjhij -\- {p \fihij F h^/j,), 

— hij^jhj]^ -\- hji^Qjhij F hij{p Pjhij F hji^{p l)/ij^ , 


F hijQjhjf^ F hj].^Qjhij F hijQjhij F hji^Qjhj ^^, 
= {hij F hji^^Qj{hij F hjf^^ F 0. 


By multiplying left and right by A^ g A 

Eq. (ID^. 


j yy Ak, we obtain 

□ 


If the terms hij are commuting, then Eq. (Dl) will 


clearly hold for any 0 < A < 1 as the left-hand side will 
be positive, while the right-hand side will be negative. In 
this case, we have the following: 


Corollary 2. Let H and H' be Hamiltonians defined as 
in Lemma^with the additional property that interaction 
terms hij commute. Then H' is gapped if 




(D7) 


Proof. As 77 is commuting, we can set A = 1 in Eq. (D51. 
Equation (D7) then follows by solving 0 < A' for p and 
restricting to p>l. □ 


We will now apply this result to our specific Hamilto¬ 
nian. Consider the map H^\5) from virtual to physical 
degrees of freedom on an injective region Ra. We can po¬ 
lar decompose = Qa{5)Wa{&) where Wa{5) is an 

isometry onto the image of and Qa{^) is an in¬ 

vertible positive operator that maps the image of Wa{5) 
onto itself. As the action of Qa{5) on (imgIUa(<5))-*- is 
arbitrary, we assume that Qa{5) is a square matrix that 
acts only the imgIUa((5). 

Let Aa(i5) = 7((5 )“^Qo(^)~^, where 'y{6) is the small¬ 
est eigenvalue of Q{6)~^ (which we have included to en¬ 
sure that A£j((5) F 1 as is required in Lemma [^. Con¬ 
sider the block tensor for two neighbouring in¬ 

jective regions Ra and Rb and let ha^b{b) be the pro¬ 
jection onto the orthogonal complement of the image of 
{Aa{5)®Ab{5))B^‘^''>\6). Themap A,(5)H(“)((5) oc Wa{5) 
is proportional to an isometry and therefore is propor¬ 
tional to 77(“)(0) up to left-multiplication by a unitary. 





















14 


Therefore, the map {Aa{S)^Ab{S))B^°‘'’^\S) is equivalent 
to up to a product of two unitaries, one act¬ 

ing on Ra and one acting on Rb. As the interaction terms 
n(“’^)((5) (the projectors onto the orthogonal complement 
of the image of [0)) are pairwise commuting, as 

we showed in the proof of Lemma[^ so are ha^b{S)- Thus 
the Hamiltonian J2{a b) ^a,b{5) is commuting and gapped, 
with a system-size independent gap of A = 1. 

Now consider the terms 

K b{5) = (A,((5) 0 Ab{5)) K,b{5) {Aa{5) ® A6(J)) , 

(D8) 

and the Hamiltonian h’^ ^(5). We see that ha,b{S) 

and f,((5) have the properties required of hij and in 


Corollary ([^. We also have that the kernel of ha,b{5) co¬ 
incides with the kernel of the projector (i.e., the 

image of B^°-A )Thus for each edge (a, h) we have that 
h'abi^) — ?7((5)n(“’^)(^), where r]{S) is the largest eigen¬ 
value of h'^ Therefore, H^{d) = J2{ab) is 

gapped for any d where b) ^'a bi^) gapped. From 
Corollary [^we thus conclude that H^{6) is gapped when 
the larg est e igenvalue /i(i5) of Aa{6) satisfies the inequal¬ 
ity Eq. ( |D7[ ). Note that from the definition of Aa(<5) that 
/r(<5)“^ can be computed as the smallest non-zero eigen¬ 
value of {B^^‘^\5)Y {S), which can be expressed as 

a simple tensor contraction. 
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